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Abstract 

We  analyze  the  propagation  of  a  high-power  laser  beam  through  a  lens  and  calculate  the 
optical  beam  quality  resulting  from  geometrical  aberrations  and  thermal  nonlinearities. 

We  present  a  general  ray  optics  formulation,  including  diffraction  effects,  for  propagation 
through  a  nonlinear  medium.  An  analytical  expression  for  the  beam  quality  parameter  M 2 
is  derived  that  includes  the  lowest-order  effects  of  geometrical  and  thermal  aberrations  in  a 
thin  lens.  A  ray  optics  simulation  including  thermal  effects  is  used  to  model  propagation 
through  a  multi-lens  optical  beam  expander  of  the  type  used  in  recent  long-range,  high- 
power  fiber  laser  experiments. 
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I.  Introduction 


High  energy  laser  systems  for  directed  energy  applications,  e.g.,  high -power  fiber  lasers  [1],  are 
prone  to  nonlinearities  in  the  optics  that  degrade  the  optical  beam  quality  and  result  in  a  large 
spreading  angle.  For  high-power  lasers,  the  beam  degradation  due  to  nonlinearities  in 
transmissive  optics  may  dominate  atmospheric  effects  such  as  turbulence,  among  others  [2].  The 
objective  of  this  Paper  is  to  analyze  and  numerically  simulate  the  effects  of  geometrical  and 
thermal  aberrations  in  lenses  on  the  optical  beam  quality  of  a  high-power  laser  beam. 

Recent  high-power  fiber  laser  experiments  conducted  by  the  Naval  Research  Laboratory 
(NRL)  at  the  Naval  Surface  Warfare  Center  (NSWC,  range  1.2  km)  and  the  Starfire  Optical 
Range  (SOR,  range  3.2  km)  propagated  a  number  of  multi-kW  fiber  laser  beams  through  the 
atmosphere  and  incoherently  combined  the  beams  onto  a  10  cm  radius  target.  In  both 
experiments,  each  beam  passed  through  a  multi-lens  beam  expander  configuration.  In  the 
NSWC  experiments,  3kW  (cw)  was  propagated  with  an  efficiency  >  90%  at  a  range  of  1.2  km 
[3].  It  was  observed  that,  as  the  power  of  an  individual  fiber  laser  was  increased  to  ~1  kW,  the 
laser  spot  size  on  target  expanded  and  was  deflected,  and  the  propagation  efficiency  dropped. 
Analysis  and  simulations  indicate  that,  at  these  power  levels  and  propagation  geometries,  thermal 
blooming  in  the  atmosphere  can  be  mitigated  by  using  a  fan  near  the  laser  source  to  induce  air 
flow  across  the  beam  path.  However,  because  of  high  laser  intensity  emerging  from  the  fiber, 
thermal  nonlinearities  in  the  transmissive  optics  cannot  be  mitigated  in  this  manner  and  result  in 
the  observed  beam  degradation. 

In  this  Paper,  we  analyze  the  propagation  physics  of  a  high-power  laser  beam  through  an 
optical  beam  expander  consisting  of  a  system  of  expanding  and  focusing  lenses.  In  particular,  we 
analyze  the  optical  beam  quality  resulting  from  geometrical  aberrations  and  thermal 
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nonlinearities.  Section  II  presents  the  general  ray  optics  formulation  for  propagation  through  a 
nonlinear  lens.  We  discuss  how  to  model  diffraction  effects  using  a  ray  optics  approach  and 
define  the  optical  beam  quality  parameter  M 2  [2].  In  Section  III  we  discuss  potential  sources  of 
nonlinearity  in  a  lens  and  analyze  the  lowest  order  effects  of  spherical  and  thermal  nonlinearities 
in  a  thin  lens  on  beam  quality.  Section  IV  presents  results  of  ray  optics  simulations  on  systems 
of  lenses.  We  examine  thermal  effects  in  a  single,  optically  thick  lens  and  model  a  multi-lens 
beam  expander  similar  to  that  used  in  recent  NRL  experiments.  In  Section  V,  we  discuss  our 
theoretical  findings  in  the  context  of  recent  experimental  results  and  present  conclusions  in 
Section  VI. 


II.  Ray  Optics  Formulation 

The  ray  optics  equations  describing  the  propagation  of  a  ray  through  a  medium  with  a  non- 
uniform  refractive  index  n{  r)  are  dr/dt  =  dco/dk,  and  dk/dt  =  -d  col  dr,  where  ris  the 

transverse  position  of  the  ray  and  k  =  |k|  is  the  magnitude  of  the  wave  vector.  The  angular 
radiation  frequency  co  =  ck/n( r)  serves  as  the  Hamiltonian  and  dco!dt  =  0  [4], 


i)  Ray  Equations 


Transforming  to  the  independent  variable  z  ,  the  ray  optics  equations  can  be  combined  to  give 
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where  x(z,t )  and  y(z,t)  denote  the  transverse  coordinates  of  the  ray  and  n(x(z),y(z),z,t)  is 
the  refractive  index  evaluated  at  the  position  of  the  ray  at  time  t .  Thermal  effects  cause  the 
refractive  index  to  vary  with  time.  In  deriving  Eqs.  (1),  we  have  made  the  quasi-static 
approximation,  i.e.,  it  is  assumed  that  the  index  varies  on  a  time  scale  long  compared  with  the 
transit  time  of  a  photon  through  the  medium.  Hence,  terms  proportional  to  the  time  derivative  of 
the  refractive  index  are  neglected,  but  the  parametric  time  dependence  of  the  refractive  index  is 
retained. 


For  axial  symmetry  with  respect  to  the  z-axis,  Eqs.  (1)  can  be  combined  to  give 


d2r_l  ^Br''  dn  dr  dii^ 
dz2  n  J^3r  dzdzj’ 


(2) 


where  r  (z)  is  the  radial  position  of  a  ray  at  axial  location  z  . 


ii)  Ray  Distribution 

The  propagation  of  a  laser  beam  through  a  refractive  medium  can  be  described  by  a  distribution 
of  N  rays  in  the  phase  space  (r,  v) ,  i.e., 

N 

/( r,  v,  z)  =£  <?(r  -  r  (r0i ,  v0i ,  z))S(\  -  v(r0i ,  v0|. ,  z))  (3) 

;=i 

where  r  and  v  denote  transverse  ray  position  and  transverse  velocity,  respectively.  The 
quantities  r(r0;,  v0i, z)  and  v(r0l , \0i , z)  are  the  phase  space  trajectories  of  the  ith  ray  with  phase 

space  coordinates  (r0i,  v0()  at  z  =  0.  These  trajectories  evolve  according  to  Eqs.  (1)  or  (2).  In 
the  limit  of  a  large  number  of  rays,  a  continuous  distribution  can  be  formed,  i.e., 

/( r,  v, z)  =  jd2r0d2\0  /( r0 ,  v0 ,0)  <?(r -  r (r0 ,  v0 ,  z))5{ v  -  v(r0 ,  v0 , z)) .  (4) 
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The  number  density  of  rays,  Jd2v/(r,  v,z) ,  is  proportional  to  the  intensity.  For  an  axially 
symmetric  laser  beam,  the  intensity  at  position  (r,  z)  is  given  by 

I(r,z)=  jl(r0,0)g(v0)—  r^°,vQ’^  r0dr0d2\0,  (5) 

where  the  subscript  zero  denotes  an  initial  (z  =  0)  value,  and  g(v0)  is  the  initial  transverse 
velocity  distribution  normalized  so  that  Jg(v0)d2v0  =  1 . 


iii)  Optical  Beam  Quality 

/  7  \1/2 

The  laser  spot  size  is  defined  as  R(z)  =  {7  (r0,\0,z))  ,  where  the  averaging  operator  for  any 

quantity  Q  is  (Q)  =  (4n/P)jQ(r0, \0,z)l(r0,0)g(\0)r0dr0d2v0 , 

and  P  is  the  beam  power.  It  can  be  shown  that  the  spot  size  satisfies  /?’  \z)  =  £2 /R3(z),  where 

e2  =(r2)((r')2  +rr")-(rr')2  (6) 

and  de2  /dz  =(r  2)(3r'  7 "  +  77”).  The  quantity  £  is  the  optical  equivalent  of  particle  beam 

emittance  [5].  In  a  uniform  medium,  r"=0  and  £2  is  a  constant  of  motion  (d£2 /dz  =  0) .  The 
optical  beam  quality  parameter  is  defined  as 

M2  =  7t£l  X  (7) 


and  is  proportional  to  the  rms  area  of  the  ray  distribution  in  optical  phase  space, 
fundamental  Gaussian  beam,  £-  XI  n  and  M 2  =  1 . 

In  a  uniform  medium,  the  equation  for  the  spot  size  has  the  solution 

f  \2 
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where  Lf  is  a  constant  that  determines  the  focal  length  and  the  asymptotic  spreading 
(diffraction)  angle  is  given  by  0  =  6/  Ro=M20o,  where  0O  =  A/(nR0)  is  the  diffraction  angle  of 

a  Gaussian  beam.  Hence,  M2  is  the  ratio  of  the  spreading  angle  of  a  beam  to  the  diffraction 
angle  of  a  Gaussian  beam  with  an  equivalent  spot  size.  The  asymptotic  spreading  angle  of  the 
optical  beam  is  referred  to  as  M 2  “times  diffraction-limited.” 

iv)  Diffraction 

The  ray  equations  are  valid  in  the  limit  where  A  — >  0  and  hence,  do  not  describe  diffraction. 
However,  diffraction  can  be  introduced  in  an  ad  hoc  fashion  by  giving  the  rays  a  transverse 
velocity  distribution  that  results  in  the  intensity  profile  remaining  Gaussian  and  having  the 
correct  vacuum  diffraction  angle. 

To  obtain  the  correct  velocity  distribution,  we  evaluate  £  for  a  Gaussian  beam  with  an 
initial  intensity  profile  7(r,z  =  0)=(2/>//rf?o)exp(-2r2//?o).  The  velocity  distribution  of  the 

rays  is  taken  to  have  the  form  g(v0)  =  (l//rv2)exp(-|  v0 12  /v2),  where  v  is  the  transverse 

velocity  spread.  A  Gaussian  velocity  distribution  ensures  that  the  intensity  profile  remains 
Gaussian  when  propagating  in  vacuum.  For  a  focused  Gaussian  beam  propagating  in  vacuum, 
Eq.  (2)  can  be  integrated  to  obtain  r  (z)  =  r0  +  [(v0  / c)-  (r0  /  Lf  )Jz .  Using  this  result  to  evaluate 

fin  Eq.  (6)  and  setting  M 2  =1  yields  the  velocity  spread  v  =  0oc/yf 2  =  Ac /(y[2nR0) ,  which 
gives  a  spreading  angle  equal  to  the  Gaussian  diffraction  angle. 

III.  Effects  of  Nonlinearities  on  Optical  Beam  Quality 

We  derive  a  general  expression  for  M 2  due  to  the  lowest  order  (dominant)  geometrical,  thermal, 
or  other  nonlinear  aberrations  in  a  thin  lens.  The  thin  lens  approximation  states  that  the  change 
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in  the  radial  position  of  a  ray  as  it  passes  through  a  lens  is  negligible.  It  is  also  assumed  that  (1) 
the  beam  is  Gaussian,  (2)  thermal  and  other  nonlinear  effects  are  proportional  to  the  beam 
intensity,  and  (3)  the  transverse  extent  of  the  nonlinear  medium  is  much  larger  than  the  laser  spot 
size.  With  these  assumptions,  the  angle  of  a  ray  after  it  passes  though  a  thin  lens  can  be  written 
as 


r'  =  ^+^-[L(r)(n(r)-l)]r=? 
or 


7  ^ 
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where  /?0  =  v0  /c  is  the  angle  of  the  ray  before  it  enters  the  lens,  L(r)  is  the  lens  thickness  and 
the  constants  a )  are  functions  of  the  geometrical  and  nonlinear  aberrations.  Using  Eq.  (9)  to 

evaluate  M2  in  Eq.  (7),  and  noting  that  (?n^  -  2l-<n/2)r(l  +  (n/2))/?o  ,  results  in 


M2  = 


n2Rl(r 

1  +  ^ 


9*' 
32  J 


2  3  \n  1 

a ;  +-J—a7a ,  +-a 
2  8  V  2  2  3  2  • 


1/2 


(10) 


The  quantities  a2  and  a3  contain  terms  describing  the  lowest  order  geometrical  aberrations  and 

other  terms  which  are  proportional  to  nonlinear  refractive  index  of  the  lens.  The  nonlinear 
refractive  index  of  a  lens  contains  contributions  from  many  physical  processes,  for  example, 
thermal  effects,  electronic  vibration  (optical  Kerr  effect),  and  electrostriction,  among  others  [6]. 

In  the  steady  state,  the  effective  nonlinear  index  due  to  thermal  effects  is  approximately 
given  by 


n2T  =  y(2R2L/R20) 
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where  RL  is  the  lens  radius,  R0  is  the  laser  spot  size,  dn/dT  is  the  thermo-optic  constant,  a  is 
the  absorption  coefficient,  k  is  the  thermal  conductivity  of  the  lens,  and  y(x)  is  a  function  of 
order  unity  that  is  defined  in  Sec.  IV.  For  a  typical  BK7  lens  with  dti/dT  =  \0~5  K~' , 
a  =  5x10^* cm-1 ,  RL  =  2.5  cm,  and  a  laser  beam  with  R0  ~  1  cm,  the  effective  nonlinear  index 
due  to  thermal  effects  is  n2T  ~  3x1  CT8  cm2/W ,  which  is  much  larger  than  the  electronic 
nonlinear  index,  i.e.,  the  optical  Kerr  effect,  for  which  n2K  ~  10“16  cm2AV  .  Hence,  from  here 
on,  we  retain  only  the  dominant  thermal  nonlinearity. 

IV.  Thermal  Nonlinearity  in  Lenses 

A  laser  beam  with  spot  size  R0  passing  through  a  lens  with  radius  RL  >  R0  will  non-uniformly 

heat  the  lens  and  modify  its  optical  properties.  We  calculate  the  steady  state  temperature  profile 
and,  in  the  following  Section,  calculate  its  effect  on  beam  quality.  We  assume  that  the 
absorption  coefficient  of  the  lens  is  small  enough  that  we  can  neglect  attenuation  of  the  laser 
within  the  lens.  The  temperature  within  the  lens  satisfies  the  heat  equation, 

p0CpdT/dt  =  rtf2T  +  al(r),  (11) 

where  p0  is  the  lens  mass  density,  Cp  is  the  specific  heat,  k  is  the  thermal  conductivity,  a  is 
the  absorption  coefficient,  and  /(r)is  the  laser  intensity.  Within  the  lens,  it  is  assumed  that  the 
heat  flux  in  the  longitudinal  direction  is  much  smaller  than  in  the  transverse  direction,  i.e.,  the 
temperature  is  uniform  in  the  longitudinal  direction.  Hence,  we  neglect  longitudinal  derivatives 
in  the  heat  equation  and  treat  temperature  as  solely  a  function  of  radial  coordinate  and  time.  In 
solving  Eq.  (11),  the  boundary  conditions  are  taken  to  be  1)  the  radial  derivative  of  the 
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temperature  is  zero  at  r  =  0 ,  and  2)  the  edge  of  the  lens  is  maintained  at  a  specified  temperature, 
T(Rl)  =  Tx,  where  RL  is  the  radius  of  the  lens. 

For  a  cw  laser  beam,  the  temperature  profile  will  approach  a  steady  state  on  a  time  scale 
rT  ~  p0CpRq/  k,  where  R0  is  the  spot  size  of  the  laser  beam.  Assuming  a  Gaussian  laser 

intensity  profile,  i.e.,  /(r)  =  /0exp(-2r2/i?02),  Eq.  (11)  has  an  analytic  steady  state  solution 
given  by 


A  T(r)  =  T(r)-T„=T 


y  ’  «»2  , 


n- 1 


f_ 

2r2  ^ 

-  r 

0, - r- 

K) 

(  2  Y 

■X 

jaT,2n 

1 

r 

UoJ 

-In 


kKj 


(12) 


where  T  =  aP KAnic) ,  aT  2n  =  (-2)"  /(nn\) ,  and  T(n,x)  is  the  incomplete  Gamma  function.  The 
steady-state  temperature  change  on-axis  is  given  by  AT (0)  =  y{lR2L  / R^)t  ,  where 
y(x)  =  y0  +  r(0,x)  +  ln(jt),  and  y0  ~  0.58  is  Euler’s  constant.  Figure  2  plots  the  radial 


temperature  profile  [Eq.  (12)]  normalized  to  T  for  various  values  of  RL/  R0. 

We  consider  the  propagation  of  a  high-power  laser  beam  through  a  lens  with  a  non-zero 
absorption  coefficient  and  spherical  surfaces  as  shown  in  Fig.  1  and  calculate  the  beam  quality 
modification  due  to  thermal  effects.  It  is  assumed  that  the  beam  spot  size  is  much  smaller  than 
the  radius  of  curvature  of  the  lens  so  that  only  the  lowest  order  spherical  aberration  is  retained. 

To  lowest  order  in  the  temperature  change,  AT ,  the  angle  of  a  ray  at  radial  position  r0 
after  passing  through  a  lens  is  given  by 

V  =  &+  o-D  +  gAT(r))]r=ro ,  (13) 
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where  J30  =  v0  /c  is  the  angle  of  the  ray  before  it  enters  the  lens,  n0  is  the  linear  refractive  index, 
g  =  (dn/dT)  +  (n0 -\)aT ,  aT  is  the  thermal  expansion  coefficient,  and  Lc(r)  is  the  lens 


thickness  in  the  absence  of  thermal  effects.  In  writing  Eq.  (13),  we  have  made  the  thin-lens 
approximation,  i.e.,  the  radial  position  of  a  ray  does  not  change  as  it  passes  though  the  lens. 

For  a  spherical  lens, 


where  Rcl  and  RC1  are  the  radii  of  curvature  of  the  two  faces  as  indicated  in  Fig.  1, 
aC2  =(R20/2L0)(R-\+R-\)and  aC4  =  (Rq  /SL0)(R^]  +R^\),  and  LQ  is  the  lens  thickness  on 
axis.  The  second  line  of  Eq.  (14)  is  an  expansion  for  r  «  RCI,RC2  which  retains  the  lowest 
order  spherical  aberration. 

In  the  analysis  that  follows  we  assume  a  Gaussian  laser  beam  and  a  bi-convex  or  bi¬ 
concave  lens  for  which  RCI  =  RC2  =  Rc ,  and  Rc  >  0  (Rc  <  0)  denotes  a  bi-concave  (bi-convex) 

lens.  For  a  bi-convex  lens,  the  lens  radius,  RL ,  is  limited  by  geometry,  i.e.. 


rl<rc 


®  (4 1  Rc  l)'/2  • 


In  the  vacuum  region  after  the  lens,  r0"=0,  so  that  e2  V)2  •  Substituting  Eqs. 


(12)  and  (14)  into  Eq.  (13)  gives 
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«,  =Kn0-\)  +  gATm^-^gT  , 


(15b) 
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Equation  (10)  then  yields  (for  a2  =0),  M 2  =  [l  +  Asa  +  Ar]l/2  where 
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are  lowest  order  modifications  to  the  beam  quality  due  to  spherical  aberration  (Eq.  (16a)),  and 
lens  heating  (Eq.  (16b)).  In  writing  Eqs.  (16),  we  have  assumed  that  the  temperature  change  is 
sufficiently  small  that  we  can  neglect  terms  that  are  of  order  ( gT )2  and  that  gAT (0) «  n0  —  1 . 

Eqs.  (16)  are  valid  when  AS/,  >>^1.  When  2 R%  /(lo|/?c|)<1,  Ar  has  the  same  sign  as  Rc, 
i.e.,  the  temperature  aberration  increases  M2  when  Rc> 0  (concave  lens)  and  decreases  M2 
when  Rc  <  0  (convex  lens).  When  2R^  /(L0|ftc|)>  1 ,  the  thermal  aberration  decreases  M 2 
independent  of  the  sign  of  Rc  . 

A  more  accurate  calculation  requires  that  more  terms  be  kept  in  the  representation  of 
temperature  since  the  expansion  parameter  r02  /  R]  >  1  for  rays  near  the  edges  of  the  intensity 

distribution.  Additionally,  one  can  relax  the  requirement  that  the  laser  spot  size  be  much  less  the 
lens  radius.  In  this  case,  it  can  be  shown  that 


ASA  =  f^hJ^hl(f), 


el 


(17a) 


and 


11 


7n(f)  =  2Mn/2)  r(l  +  (n / 2)) - if  1  +  (n / 2), y 

v  J 

r(n)  is  the  Gamma  function,  and  r(n,;c)is  the  incomplete  Gamma  function.  For  a  bi-convex 
lens,  |<rj<l.  The  sum  in  Eq.  (17b)  requires  ~0(R2LIRl)  terms  to  converge.  The  functions 
hn(f )  are  positive  definite.  Letting  /  — > 0  and  retaining  only  the  first  term  in  Eq.  (17b) 
recovers  Eqs.  (16). 

When  analyzing  Eqs.  (17)  and  in  the  numerical  simulations  that  follow,  we  will  consider 
lenses  made  of  either  BK7  glass  or  optical-grade  fused  silica.  Laser  absorption  in  fused  silica  in 
the  infrared  is  primarily  due  to  the  OH'  content,  which  can  vary  by  orders  of  magnitude 
depending  on  the  type  of  fused  silica  and  how  it  is  manufactured.  Suprasil  W-l  is  a  variant  of 
fused  silica  with  an  extremely  low  OH  content  of  <  3ppm  [7]  and  an  absorption  coefficient  ~  10'6 
cm'1  at  1  |xm  [8].  Optical-grade  fused  quartz  has  an  OH  content  of  ~  400-500  ppm  and  an 
absorption  coefficient  of  ~  10‘5  cm'1,  while  standard  UV-grade  fused  silica  has  a  higher  OH 
content  of  >  1000  ppm  [9].  The  thermo-optical  properties  of  these  two  materials,  summarized  in 
Table  1,  are  taken  from  Ref.  [10]  except  for  the  absorption  coefficients  of  BK7  and  fused  silica, 
which  were  measured  in  our  laboratory.  The  fused  silica  lens  used  in  these  measurements  had  an 
OH  content  of  >  1000  ppm.  The  value  of  dn/dT  for  fused  silica  was  estimated  from  Ref.  [10]. 

Figure  3  plots  Asa  and  AT  (Eq.  (17))  as  functions  of  the  normalized  laser  beam  spot  size, 

Rq  /  Rl  ,  for  a  bi-concave  ( Rc  >  0 )  and  bi-convex  ( Rc  <  0 )  lens.  The  lens  is  assumed  to  be  BK7 
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glass  with  n0  =  1 .45 ,  |/?c|  =  30  cm ,  RL  =  3  cm ,  and  L0  =  1  cm .  The  laser  wavelength  is  1  pm 
and  the  power  is  1  kW.  It  is  seen  that  the  thermal  aberration  increases  M2for  Rc>  0  and 
decreases  M 2  for  Rc  <  0 . 

V.  Numerical  Simulation 

We  have  developed  a  ray  optics  code  based  on  solving  Eq.  (2)  for  an  axially  symmetric 
distribution  of  rays  which  pass  through  an  arbitrarily  specified  refractive  index  n(r,z,t).  The 
time  dependence  of  n  is  due  to  heating  and  is  self-consistently  obtained  by  solving  Eq.  (11) 
using  the  boundary  conditions  discussed  in  Section  IV.  When  calculating  the  temperature 

profile,  a  Gaussian  intensity  profile  is  assumed  with  a  spot  size  given  by  R(z)  =  (r2(r0,\0,z) )'/2. 

Diffraction  is  modeled  by  assuming  a  distribution  of  transverse  velocities  as  specified  in  Section 
III.  We  assume  a  cw,  Gaussian  beam  that  turns  on  at  t  =  0  and  calculate  the  beam  quality  after 
passing  through  a  collection  of  lenses.  Although  the  input  beam  is  assumed  to  be  Gaussian,  the 
model  places  no  limitations  on  M2  as  the  beam  propagates  through  the  optical  system.  We  first 
consider  the  case  of  single  bi-convex  or  bi-concave  lenses  to  compare  with  the  analysis  of  the 
previous  Section.  We  then  model  a  beam  director  configuration  that  was  used  in  a  recent  series 
of  high-power  fiber  laser  experiments  performed  by  NRL  in  which  thermal  effects  were 
observed  [3]. 

i)  Single  Lens 

We  consider  the  case  of  a  single  bi-convex  or  bi-concave  spherical  lens  with  the  same 
parameters  used  in  plotting  Fig.  3.  The  laser  wavelength  is  1  pm  and  the  spot  size  is 

R0  =  1 .25  cm .  The  steady-state  time  constant  for  the  temperature  is  tt  ~  p0C pRl  /  K  ~  5 
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minutes.  Figure  4a  plots  the  on-axis  temperature  change,  AT(0),  as  a  function  of  time.  Figure 
4b  plots  M 2  versus  time  for  a  bi-convex  (blue  curve)  and  a  bi-concave  (red  curve)  lens.  At 
t  =  0,  in  the  absence  of  heating,  M2  =2.7  for  both  the  concave  and  convex  lenses.  Since  the 
quantity  4|<r|/  -  0.2  is  less  than  unity  for  these  parameters,  M 2  decreases  (increases)  as  the 

temperature  increases  for  the  bi-convex  (bi-concave)  lens,  consistent  with  the  analysis  of  the 
previous  Section.  The  beam  quality  reaches  a  steady  state  in  ~  5  minutes.  In  the  steady-state, 
M2  ~  3.2  for  the  concave  lens  and  M2  ~  2.2  for  the  convex  lens. 

Figure  5  plots  A/ 2  in  the  steady  state  versus  laser  power  for  different  lens  thicknesses  for 
convex  and  concave  lenses.  For  the  case  of  a  bi-concave  lens,  M 2  ~  2.7  when  the  power  is  near 
zero  and  increases  with  power  and  lens  thickness.  For  a  bi-convex  lens,  M 2  ~  2.7  when  the 
power  is  near  zero  and,  for  the  two  thinner  lenses,  decreases  as  the  beam  power  is  increased  over 
the  range  plotted.  For  the  thicker  lens,  M 2  initially  decreases  and  then  increases  with  power.  A 
minimum  M 2  of  ~1.4  is  obtained  when  the  power  is  -  4  kW. 

ii)  NRL  Beam  Expander 

The  earlier  beam  expander  configuration  used  in  the  NSWC  experiment  had  five  lenses  made  of 
various  optical  materials  such  as  fused  silica  and  BK7.  In  that  experiment,  it  was  observed  that 
beam  quality  deteriorated  severely  when  the  laser  power  approached  -  1  kW.  This  was 
attributed  to  thermal  effects  in  the  optics.  A  more  recent  version  of  the  NRL  beam  expander  is  a 
3-lens  system  made  of  uv-grade  fused  silica  that  we  model  here  to  examine  the  effects  of  thermal 
aberrations. 

A  schematic  diagram  of  the  beam  expander  is  shown  in  Fig.  6  with  the  lens  parameters 
and  geometry  listed  in  Table  1.  The  first  two  plano-concave  lenses  expand  the  beam  from  a  spot 
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size  of  1 .25  cm  to  2.5  cm.  The  short  focal  lengths  of  these  first  two  lenses  introduce  a  large 
degree  of  spherical  aberration  which  is  corrected  by  the  third  lens.  For  an  input  beam  with 
M2  =1 ,  in  the  absence  of  thermal  aberrations,  the  beam  quality  is  M2  ~  4  after  the  first  lens, 
M 2  ~  20  after  the  second  lens,  and  M 2  ~  1  after  the  third  collimating  lens.  This  low-power 
limit  has  been  successfully  benchmarked  against  the  ZEMAX  optical  design  code  [11]. 

In  the  presence  of  thermal  aberrations,  M 2  increases  with  laser  power.  Figure  7  plots  the  output 
M 2  versus  laser  power  for  a  beam  expander  made  from  BK7  lenses  (dashed  curve)  and  fused 
silica  lenses  (solid  curve).  The  thermo-optical  properties  of  each  material  are  listed  in  Table  1. 
For  simplicity,  it  is  assumed  that  the  refractive  indices  of  both  materials  are  n  =  1 .45 . 
Compared  with  fused  silica,  BK7  has  higher  absorption  and  thermal  expansion  coefficients,  but  a 
lower  thermo-optic  coefficient.  The  overall  result  is  that  M2  increases  faster  with  increasing 
laser  power  for  the  BK7  configuration,  e.g.,  for  a  laser  power  of  ~  1.5  kW,  M2  ~  3  for  the  BK7 
lenses  and  M2  ~  1 .2  for  fused  silica  lenses. 

In  Figure  8,  we  start  with  a  beam  expander  with  three  fused  silica  lenses,  replace  one  of 
the  lenses  with  BK7,  and  plot  M 2  versus  laser  power.  It  is  seen  that  the  beam  quality  has  a 
stronger  dependence  on  laser  power  ( M 2  increases  faster  with  laser  power)  when  it  is  the  first 
lens  that  is  replaced  with  BK7.  The  weakest  dependence  results  from  replacing  the  third  lens. 
Comparing  with  Fig.  7,  we  note  that  when  the  first  lens  is  BK7  and  the  second  and  third  lenses 
are  fused  silica,  the  resulting  M2  is  larger  than  that  of  case  where  all  the  lenses  are  BK7.  Also, 
the  beam  expander  in  which  lenses  1  and  2  are  fused  silica  and  lens  3  is  BK7  has  a  lower  M 2 
than  a  beam  expander  in  which  all  of  the  lenses  are  made  of  fused  silica. 
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VI.  Discussion 


We  calculate  the  laser  propagation  efficiency  for  the  conditions  of  the  NRL  experiment  at  SOR 
and  compare  with  the  observations.  The  propagation  efficiency  is  defined  as  the  ratio  of  the 
power  on  a  target  to  the  transmitted  power.  For  a  single  Gaussian  beam  with  power  P0  and  a 

target  of  radius  /?target  at  range  L,  the  propagation  efficiency  is  given  by 

V  =  Pm get  IP,  =exp (-J3L)  [l  -  exp(-2/?^get  IR2(L))\,  where  Pmget  is  the  power  on  the  target  and 


[3] 


R(L)  =  R( 0) 
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(18) 


is  the  beam  spot  size  on  the  target,  R( 0)  is  the  spot  size  after  exiting  the  beam  expander. 


P„  =  0A58(X2  l(C2n  L)f5  is  the  transverse  coherence  length  due  to  turbulence,  and  C2  is  the 
turbulence  strength  parameter  [2].  In  writing  Eq.  (18),  it  is  assumed  that  the  focal  length  is  equal 
to  the  range  to  the  target  and  that  M2  ~  1 .  Equation  (18)  becomes  less  accurate  when  M2  » 1 
because  of  the  coupling  of  higher  order  modes  with  turbulence. 

The  propagation  efficiency  is  plotted  in  Fig.  9  as  a  function  of  laser  power  for  the 
conditions  of  the  experiment,  i.e.,  X  —  \pva,  C 2  =10_l5m'2/3,  L  =  3km,  /?  =  0.08km'1,  and 
^target  =  10  cm  •  It  is  seen  that  for  BK7  lenses,  the  propagation  efficiency  decreases  rapidly  with 
laser  power  when  P0  >  500  W  .  For  fused  silica,  the  propagation  efficiency  remains  relatively 

constant  (~  80%)  for  laser  powers  <  2  kW.  Hence,  a  kW-class  beam  director  using  transmissive 
optics  will  require  the  use  of  a  low-absorption  material  such  as  fused  silica. 

The  value  of  M 2  from  the  laser  is  measured  in  the  laboratory  to  be  1.1  ±0.05  at  1  kW  of 
power.  It  is  observed  that  M 2  increases  by  30%  after  passing  through  the  beam  expander  using 
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all  fused  silica  lenses.  Theoretically,  from  Fig.  7,  M 2  increases  by  ~20%  at  1  kW.  There  are 
several  differences  between  the  beam  expander  modeled  here,  and  the  actual  beam  expander 
used  in  the  experiments  that  may  account  for  this  difference.  First,  we  have  neglected  the 
heating  and  convection  of  the  air  between  the  lenses.  The  beam  expander  lenses,  when  used  in 
the  experiment,  were  sealed  inside  a  metal  mounting  tube  which  could  enhance  the  heating  of  the 
air  surrounding  the  lenses.  Second,  lens  misalignment  due  to  heating  of  the  lens  mounts  is  not 
accounted  for  in  the  model.  Third,  the  theoretical  radial  boundary  condition  on  the  temperature 
profile  may  not  be  accurate  for  laser  spot  sizes  comparable  to  the  lens  radius.  Using  a  boundary 
condition  that  assumes  radiative  loss  at  the  edge  of  the  lens  results  in  a  larger  temperature 
change.  However,  in  this  parameter  regime,  the  beam  quality  is  not  affected  significantly  by  the 
choice  of  boundary  conditions.  Fourth,  there  is  uncertainty  in  our  experimental  measurements  of 
the  absorption  coefficient  of  the  lenses.  Lastly,  we  note  that  significant  beam  degradation  occurs 
within  the  laser  source  itself  when  it  is  run  at  higher  power  (~  2  kW). 

VII.  Conclusions 

We  have  analyzed  and  modeled  the  propagation  of  a  high-power  laser  beam  through  a  lens  and 
calculated  the  optical  beam  quality  resulting  from  geometrical  aberrations  and  thermal 
nonlinearities.  It  is  found  that  for  a  thin,  spherical  lens  in  the  low-power  regime,  thermal 
aberrations  can  either  enhance  or  degrade  the  beam  quality  depending  on  the  value  of  the 
parameter  a  =  2 R%  /(lo|/?c|)  and  the  curvature  of  the  lens  (convex  or  concave)  .  When  a  <  1 , 

the  temperature  aberration  increases  M 2  for  concave  lens  and  decreases  M 2  for  a  convex  lens. 
When  a  >  1 ,  the  thermal  aberration  decreases  M 2  independent  of  the  sign  of  Rc . 
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We  have  modeled  a  3-lens  optical  beam  expander  similar  to  that  used  in  recent  high- 
power  fiber  laser  experiments  performed  by  NRL.  It  is  found  that  when  using  kilowatt-class 
lasers  and  transmissive  optics,  low-absorption  material  such  as  optical-grade  fused  silica  is 
required  to  achieve  high  propagation  efficiencies  over  kilometer  ranges. 
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Material 

BK7 

Fused  Silica 

a :  Absorption  coefficient 

3x10^  cm'1 

5xl0"5  cm-1 

K\  Thermal  conductivity 

l.lxlO-2  W/(cm K) 

1 .4  x  1 0“2  W/(cm  K) 

dn/dT 

2.3X10-6  K~' 

o 

1 

ar=L~ldL/dT 

7-lxlO-6  K ~x 

5.1xl0“7  K~' 

p0 :  density 

2.5  g/cm3 

2.2  g/cm3 

Cp\  heat  capacity 

0.86  J/(g  K) 

0.75  J/(g  K) 

Table  1:  Thermal  and  optical  parameters  for  BK7  and  uv-grade  fused  silica  used  in  the  models. 
Parameters  are  taken  from  Ref.  [10],  except  for  the  absorption  coefficients  of  BK7  and  fused 
silica,  which  were  measured  in  the  lab.  The  value  of  dn/dT  for  fused  silica  was  estimated  from 
Ref.  [10]. 


20 


Lens  1 

Lens  2 

Lens  3 

Position,  h 

0 

1.5856  cm 

17.59  cm 

Thickness,  L, 

0.62  cm 

0.62  cm 

2  cm 

Curvature 

RCi  =18.03  cm 

RC2  =15.45  cm 

/?C31  =  -64.607  cm 

RC32  =  -20.7975  cm 

Table  2:  Lens  geometry  parameters  for  the  3-lens  NRL  beam  expander.  Position,  zt ,  denotes  the 
z-coordinate  of  the  left  face  of  the  lens.  Lt  denotes  the  on-axis  thickness.  For  lens  3,  RC3X 
( RC32 )  is  the  radius  of  curvature  of  the  left  (right)  face.  The  radial  dependence  of  the  lens 
thickness  in  the  absence  of  thermal  effects  has  the  form  of  Eq.  (14). 
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Figure  1:  Schematic  diagram  of  (a)  bi-concave  and  (b)  bi-con  vex  spherical  lenses  with  linear 
refractive  index  nL(r,z )  and  radii  of  curvature  Rcl  and  RC2  The  quantity  AL(r)  denotes  the 

lens  thickness  at  radial  position  r .  A  laser  beam  (ray  distribution)  with  an  initial  spot  size  R0 , 
transverse  velocity  spread  v ,  and  M 2  =  1  propagates  along  the  z  axis. 
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Figure  2:  Normalized  steady  state  temperature  change  [Eq.  (10)]  versus  normalized  radial 
coordinate  for  RL / R0  =  2  (red),  4  (green),  and  8  (blue),  where  RL  is  the  lens  radius,  R0  is  the 

laser  spot  size,  and  AT  =  T(r)-Tm. 
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Figure  3:  Beam  quality  modification  due  to  spherical  aberration  (AS/,)  and  thermal  aberration 
(  At  )  versus  laser  beam  spot  size.  The  laser  spot  size,  R0 ,  is  normalized  to  the  lens  radius,  RL  . 
The  beam  quality  is  given  by  M2  =  [l  +  Asa  +  Ar]'/2.  Curves  denote  Asa  (black  curve,  Eq. 
(13))  and  10  Ar  (Eq.  (14))  for  a  bi-concave  (red  curve)  and  a  bi-convex  (blue  curve)  lens.  A 
BK7  lens  is  assumed  with  n0  =  1 .45 ,  |/?c|  =  30  cm ,  RL  =  3  cm  ,  L0  =  1  cm .  The  laser 
wavelength  is  1  pm  and  the  power  is  1  kW. 
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time  [min] 


Figure  4:  Peak  temperature  (Fig.  a)  and  beam  quality  (Fig.  b)  versus  time  for  a  bi-concave  (red 
curve)  and  bi-con  vex  (blue  curve)  lens.  The  lens  parameters  are  the  same  as  in  Fig.  3.  The  laser 
wavelength  is  1  pm ,  R0  =  1 .25  cm ,  and  the  power  is  2  kW. 
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_  b)  bi-convex  lens 
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Figure  5:  Steady  state  beam  quality,  M 2 ,  versus  laser  beam  power  for  a  single  (a)  bi-concave 
and  (b)  bi-convex  BK7  lens  with  thickness  L0  =  0.5  cm  (red  curves),  L0  =  1  cm  (blue  curves), 

and  L0  =  2  cm  (black  curves).  The  other  lens  parameters  are  the  same  as  in  Fig.  3. 
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Z\  Z2  Z.  3 


Figure  6:  Schematic  diagram  of  the  NRL  high-power  fiber  laser  beam  expander.  Lens  geometry 
is  listed  in  Table  2.  A  kW-class  beam  of  diameter  ~  2.5  cm  is  expanded  to  a  diameter  of  ~  4.8 
for  propagation  in  the  atmosphere.  Lens  parameters  are  listed  in  Table  1 .  For  an  input  beam  wit 
M 2  ~  1 ,  the  beam  quality  is  M 2  ~  4  after  the  first  lens,  M 2  ~  20  after  the  second  lens,  and 
M 2  ~  1  after  the  third  lens,  in  the  absence  of  thermal  effects. 
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Figure  7:  (a)  Steady  state  beam  quality  versus  laser  power  at  the  output  of  the  expander  shown 
in  Fig.  6.  Fused  silica  (FS)  and  BK7  lenses  are  assumed  with  parameters  given  in  Table  1. 
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Figure  8:  Steady-state  M 2  versus  power  for  the  NRL  beam  expander  with  one  BK7  lens  and 
two  fused  silica  lenses.  Lens  parameters  are  listed  in  Table  1.  Curves  denote  cases  where  lens  1 
(short  dashes),  2  (longer  dashes),  or  3  (solid  curve)  is  BK7.  The  laser  wavelength  is  1  pm,  and 
the  initial  spot  size  R0  =  1 .25  cm . 
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Laser  Power  [kW] 

Figure  9:  Propagation  efficiency  versus  laser  power  for  the  beam  expander  configuration  of  Fig. 
6  with  BK7  lenses  (dashed  curve)  and  fused  silica  lenses  (solid  curves).  Propagation  parameters 
are  typical  of  the  NRL  SOR  experiments,  i.e.,  X  =  l//m ,  C2  =  10“15  m_2/3 , 

L  =  3.2  km  ,  p  =  0.08km'1 ,  and  target  radius,  RMget  =  10  cm . 
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